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A relativistic light front formulation of nuclear dynamics is developed and 
applied to treating infinite nuclear matter in a method which includes the cor- 
relations of pairs of nucleons: this is light front Brueckner theory. We start 
with a hadronic meson-baryon Lagrangian that is consistent with chiral sym- 
metry. This is used to obtain a light front version of a one-boson-exchange 
nucleon-nucleon potential (OBEP). The accuracy of our description of the 
nucleon-nucleon (NN) data is good, and similar to that of other relativistic 
OBEP models. We derive, within the light front formalism, the Hartree-Fock 
and Brueckner Hartree-Fock equations. Applying our light front OBEP, the 
nuclear matter saturation properties are reasonably well reproduced. We ob- 
tain a value of the compressibility, 180 MeV, that is smaller than that of 
alternative relativistic approaches to nuclear matter in which the compress- 
ibility usually comes out too large. Because the derivation starts from a 
meson-baryon Lagrangian, we are able to show that replacing the meson de- 
grees of freedom by a NN interaction is a consistent approximation, and the 
formalism allows one to calculate corrections to this approximation in a well- 
organized manner. The simplicity of the vacuum in our light front approach 
is an important feature in allowing the derivations to proceed. The mesonic 
Fock space components of the nuclear wave function are obtained also, and 
aspects of the meson and nucleon plus-momentum distribution functions are 
computed. We find that there are about 0.05 excess pions per nucleon. 
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I. INTRODUCTION 



We introduce a light front formalism for infinite nuclear matter, in which the effects 
of correlations are taken into account. This is a light front Brueckner theory of nuclear 
matter. While the ultimate goal of this and related studies is to provide a fully relativistic 
treatment of nuclei which includes all previous knowledge about nuclear dynamics, the 
present work represents one step beyond the previous light-front mean-field calculation |l],|2j 
of the properties of infinite nuclear matter. 

Understanding an important class of experiments seems to require that light-front dy- 
namics and the related light cone variables be used. Consider the lepton-nucleus deep 
inelastic scattering experiments || which showed that there is a significant difference be- 
tween the parton distributions of free nucleons and nucleons in a nucleus. This difference 
can interpreted as a small (~ 10%) shift in the momentum distribution of valence quarks 
towards smaller values of the Bjorken variable Xsj- The Bjorken variable is a ratio of the 
plus-momentum k + = k° + k 3 of a quark to that of the target. If one regards the nucleus as 
a collection of nucleons, Xbj = p + /k + , where k + is the plus momentum of a nucleon bound 
in the nucleus. If one uses k° + k 3 as a momentum variable the corresponding canonical 
spatial variable 3 and the time variable is x° + x 3 Q . This is the light front 

(LF) approach of Dirac |J; see the recent reviews ||[7]] for more information. 

Deep inelastic scattering depends on the light-front momentum distribution which is the 
probability f{k + ) that a bound nucleon has a momentum k + . Other nuclear reactions, such 
as (e,e') and (p,2p) depend also on this very same probability PHlO)]- The quantity f{k + ) is 
simply related to the square of the ground state wave function, computed using light front 
dynamics. The usual equal time approach to nuclear dynamics is very successful, and it is 
natural to use this information to calculate the distribution f(k + ). However, in the standard 
equal time formulation this quantity is a response function and depends on matrix elements 
between the ground and all excited states, and therefore can be more difficult to compute. 

The use of light front variables is convenient for interpreting certain experiments, but 
does not allow one to avoid the necessary task of handling nuclear dynamics. Thus one is 
faced with the task of computing the ground state nuclear wave function using x + = x° + x 3 
as a time variable. The present effort is a simplification in that the nucleus is taken to 
be infinite nuclear matter. However, the detailed effects of the interactions between two 
nucleons are included, so that we are concerned with the relativistic dynamics of a strongly 
interacting many body system. 

Light front techniques have previously been applied to systems of two hadrons 
Pj8|j9|,pl]-[r5| . Our emphasis here is in large nuclear systems. The light front quantiza- 
tion procedure necessary to treat nucleon interactions with scalar and vector mesons was 
derived by Soper JT(| , and by Yan and collaborators |T7|JT8[ . 



We next outline our procedure. The necessary Lagrangian, which respects chiral symme- 
try, and its light front Hamiltonian is described in Sect. II. Its application to nucleon-nucleon 
scattering in the one-boson exchange approximation is carried out in Sect. III. A new fea- 
ture is that the effects of isovector mesons such as the p and d, and the p-nucleon tensor 
interaction are included. The NN potential is generated using the one-boson exchange ap- 
proximation. The Weinberg-type integral equation, which maintains unitarity and boost 
invariance in one direction, is solved and the results are compared with phase shift data. 
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Sect. IV is concerned with the many-nucleon problem. Two separate perturbation series 
are involved. The first step is to eliminate temporarily the meson degrees of freedom in 
favor of our two nucleon potential. Thus one first proceeds using only nucleon degrees of 
freedom. The light-front formalism to obtain the nucleonic interacting ground state wave 
function | $) in terms of a series in which a Brueckner G matrix acts on a best Slater 
determinant | <fi) is developed. The independent-pair approximation is used. One must 
find an eigenstate of the P~ operator for which the expectation value of P + is equal to 
the eigenvalue of P~. This formalism is applied and the results are discussed in Sect. V. 
The full nuclear wave function | \^), including the meson degrees of freedom, is discussed 
in Sect. VI. The object | is related to | <3?) by a second series involving the difference 
between the nucleon-nucleon interaction and the meson nucleon interactions. One finds that 
the expression for the nuclear mass, evaluated in Sect. V, is valid within our approximation. 
Furthermore, expressions for the meson and nucleon distribution functions are obtained. We 
derive a sum rule for obtaining the total number of (non-vector) mesons in the nucleus. A 
brief discussion of the implications of our results for lepton-nucleus deep inelastic scattering 
and the nuclear Drell-Yan process is presented in Sect. VII. A brief summary of our results is 
contained in Set. VIII. Some of the necessary notation is discussed in an Appendix. Some of 
the present results, but none of the details of the derivation or of our two-nucleon potential 
have appeared in Ref. [|19|. 



II. LIGHT FRONT QUANTIZATION: LAGRANGIAN, FIELD EQUATIONS, 

AND LIGHT-FRONT HAMILTONIAN 

The light front approach is a three-dimensional formalism involving a Hamiltonian which 
is a P~ operator. One starts with a Lagrangian and derives field equations which allow one to 
eliminate the appearance of dependent degrees of freedom in the Hamiltonian. Our starting 
point is a non-linear chiral model in which the nuclear constituents are nucleons ip (or ip'), 



pions 7r, scalar mesons [20], and vector mesons V^. The Lagrangian C is given by 



C = \{d^ - mltf) - \V^V, V + + \fTr{d, U &> 17+) 

+ \m 2 J 2 Tr(U + £/t - 2) + $ (j^- d, -g v V,) - U{M + grff) $ (2.1) 

where the bare masses of the nucleon, scalar and vector mesons are given by M, m s , m v , and 
ynu _ Qnyv _ Qvyn _ ^he un jtary matrix U can be chosen from amongst three forms Iff. 

U x =j»™lt, ^ * + *75T • 7T/2/ U3 = ^_ n 2 /f 2 + , 7sT . (2 . 2) 

1 - 27 5 T • 71-/2/ v 

which correspond to different definitions of the fields. This Lagrangian, based on the lin- 
ear representations of chiral symmetry used by Gursey fl2l| , is discussed in Ref. It is 



approximately (m^ ^ 0) invariant under the chiral transformation 

i// -> e* 75X 'V, U -> e -^ 5T - a U e - il5T ' a . (2.3) 
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This invariance shows that our scalar meson is not a chiral partner of the pion. Note the 
presence of the term U(M + g s <f>) which was incorrectly given as MU + g s (j) in Refs. 

The constant M- plays the role of the bare pion-nucleon coupling constant. If / is chosen 
to be the pion decay constant, the Goldberger-Trieman relation says that the axial vector 
coupling constant g& — 1. This is not really a problem because loop effects can make up 
the needed 25% effect. Corrections of that size are typical of order (y-) 3 effects found in the 
cloudy bag model for many observables, including g&. We also note that the A is not 
treated as an explicit degree of freedom in the above Lagrangian. 

The present Lagrangian may be thought of as a low energy effective theory for nuclei 
under normal conditions. A more sophisticated Lagrangian is reviewed in ]23[ and used in 
the present one is used to show that light front techniques can be applied to hadronic 



theories relevant for nuclear physics. This hadronic model, when evaluated in mean field 
approximation, gives |25| at least a qualitatively good description of many (but not all) 



nuclear properties and reactions. There are a variety of problems occurring when higher 
order terms are included [p3 j. The aim here is to use a reasonably sophisticated Lagrangian 
to study the effects that one might obtain by using a light front formulation. 

Ref. contains the details of the quantization procedure; we re-state the relevant results 
here. An essential feature is the quantization of spin 1/2 fermions. Although described by 
four- component spinors, these fields have only two independent degrees of freedom. The 
light front formalism allows a convenient separation of dependent and independent variables 
via the projection operators A± = 7°7 ± /2 f[6 |, with ip± = A±ip±. The independent Fermion 



degree of freedom is chosen to be ip' + . The properties of the projection operators are discussed 
in Appendix A. One gets coupled equations for ij}' ± : 

(id- - g v V~W + = (a ± ■ (p ± - g v V ± ) + (3U(M + g s <p))tfj'_ 
[id + - g v V + W_ = (ol • ( P± - g v V ± ) + (3U(M + g s <t>)W + . (2.4) 

The relation between and is very complicated unless one may set the plus component 
of the vector field to zero 0. This is immediately obtained in QED and QCD by choosing 
an appropriate gauge in which the plus-component of the vector potential vanishes. Here 
the non-zero mass of the vector meson prevents such a choice. Instead, one simplifies the 



equation for ip'_ by [16,18| transforming the Fermion field according to 

if>' = e _ift,A(x V (2.5) 
with d + A = V + . This transformation leads to the result 

{id- - g v V-)^ + = (ol • (p± - g v V ± ) + (3U(M + g s <P))^. 

id + i>- = (a ± ■ ( P± - g v V ± ) + /3U(M + g s 4>))^+ (2.6) 



where 



d+V» = d + V» - d»V + = V + ^. (2.7) 



The fields V 1 enter the meson field equations, but the fields V 1 enter the fermion field 
equations. The eigenmode expansion for V^ 1 is given by: 



4 



V{x) 



d 2 k ± dk+6(k+) _ uf1 . 



where the polarization vectors e fl (k, 00) are given by 



w) = e"(k,cu) - —e + (k,uj) 
k + 



with the properties 



777- 

k%(k,co) = -jfe + (k,uj), £ ^(k,tuy(k,tu) = -(g^-g 



oj=1,3 



k V 

k+ 
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k+- 



(2.8) 



(2.9) 



(2.10) 



The use of the Fermion field equation allows one to obtain the light front Hamiltonian 
density 



V ±( j> + my + -(V- 
2 



)2 + _yklykl + m 2ykyk 



+ (V ± 7T) 2 + 



7T 



7T- 1 



1 " ~ Mn 7. 



7T 



7T 



m^/ sin - 



+2^ + [i-d -g v V~)^ + . 



(2.11) 



The expression ( [2.11|) is useful for situations, such as in the mean field approximation, for 
which a simple expression for ip + is known. This is not always the case, so it is worthwhile 
to use the Dirac equation to express T +_ in an alternate form: 



T+- = V ± <f) ■ V ± <f) + ml<j) 2 + ^(V r+ ~) 2 + W kl V kl + m 2 v V k V 



+ (V ±7 r) 2 + i^-, 

7T Z 



7T 



./ 



,7T 



1 J Z " \ 1 2 £2 • 2 

1 7 + m 7r/ S111 



7T 



7T 
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(2.12) 
. It is 



+V> (7± • (P± - 9vV±) + U(M + g 8< j>)) i). 

The relationship between ip and C, contains no interactions) is discussed in Ref. 
£ that is expanded in creation and destruction operators according to 

= / d [^2°J^ £ H k > ^' ik ' X Kk, A) + v(k, X)e^(k, A)] . (2.13) 

The spinors u(k, A) are the usual equal time Dirac spinors, of normalization uu = 2M. 
It is legitimate to use these because one is free to choose the representation of the solutions 
of the Dirac equation in an infinite number of ways. In particular, the correct fermionic 
anti-commutation relation for £ + is obtained with these spinors [lTj . 

The Hamiltonian is a sum of a free Pq(N) and interacting terms Pf(N): 



Po(N) = \j d 2 x L dx-Ul±. • p + M) £. 



(2.14) 
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with 



Pj = V i + v 2 + v 3 , 
vi = I d 2 x ± dx-£ (g vl ■ V + M(U - 1) + g a <f>u) £, 



(2.15) 
(2.16) 



v 2 = J d 2 x L dx~l (-g vl ■ V + M(U — 1) + g s <t>U) ^ (-g v i • V + M(U - 1) + g 8 <j>u) f , 

(2.17) 

and 

Vs = 32/ d2x±dx ~ J d yi£( x ±iyih + £( x ±iyi) e ( x ~ -vT) J d V2 e ( x ~ -y^)^(a ; ±,^)7 + ^(a ; ±, 

(2.18) 

The term v\ accounts for the emission or absorption of a single vector or scalar meson, as 
well as the emission or absorption of any number of pions through the operator U — 1. The 
term v 2 includes contact terms in which there is propagation of an instantaneous fermion. 
The term V3 accounts for the propagation of an instantaneous vector meson. 

The component that is related to the plus momentum is T ++ . The necessary expression 
is given by 

T++ = v ik V ik + m 2 v V + V + + $ 7 + zffV 

+<9 + 0<9 + + d + TZ ■ d + 7Z + 77 • fl + 7T , ^ (1 - ^Sm 2 ^). (2.19) 



III. NUCLEON-NUCLEON SCATTERING VIA ONE-BOSON EXCHANGE 

POTENTIALS 

The correlations between nucleons are caused by the nucleon-nucleon interaction. Thus 
a necessary first step towards a light-front theory of nuclear correlations is the derivation 
of a light-front theory of the nucleon-nucleon interaction. Previous work showed that 
the light-front version of the Lippmann-Schwinger equation, the Weinberg equation, can 
be transformed (with one difference remaining [2f|) into the Blankenbecler-Sugar equation. 



Kinematic invariance under boosts in the three-direction is maintained, and we shall obtain 
a one-boson exchange potential which is in reasonably good agreement with the NN phase 
shifts. 



A. General Formalism 

It is worthwhile to begin by reviewing how using the light-front Hamiltonian of Eqs. ( [2.14 - 
|2.18| ) leads to the one-boson exchange potential. This derivation is useful in understanding 
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how the full nuclear wave function discussed in Sect. VI is related to the nucleonic truncation 
of Sect. IV. Consider the scattering process 1 + 2^3 + 4. The use of second-order per- 
turbation theory shows that the lowest-order contribution to the nucleon-nucleon scattering 
amplitude is given by 



(3, 4\K\1, 2) = (3, ^(i^M + v 3 \l, 2), (3.1) 

with 



QoiPij) = 75—7^ , (3-2) 

where P~j is the negative component of the total initial or final momentum which are the 
same. In constructing the NN potential, one uses conservation of four-momentum between 
the initial and final NN states. The expression (|3.1|) yields a one-boson exchange approxi- 
mation to the nucleon-nucleon potential. 

It is worthwhile to discuss the energy denominator P^ — Pq in more detail. To be 
specific, suppose that kf > k 3 . Then the emitted meson of mass /x has momentum k with 
k + = kf — &3~, k± = ki ± — k 3± and k~ = fcj , + M . Then 



k+ 

k± 2 + fj, 2 

h + — h~ 



Pij Pq — A 2 Pq — Pu Pq — (^i ^3 ) ,„+ /„+ ' (3-3) 



The interaction K also contains a factor of k + in the denominator, so that the relevant 
denominator is D = k + (P^ — P ) = (k+ — k 3 )(ki — k 3 ) — k\ — fi 2 = q 2 — /i 2 . This last 
familiar form involves the four-momentum transfer between nucleons 1 and 3 (kf > kj) 
and leads to the usual Yukawa-type potentials. It is also useful to explore the form of the 
energy denominator using light-front variables by first defining the plus-component, P + , of 
the initial and final total momentum. We may also define kf = xP + and k 3 = x'P + in 
which x, x' are invariant under Lorentz transformations in the three-direction. Then using 
(p. 3D, we find 



x x 

This quantity is also invariant under Lorentz transformations in the three-direction. 
A straightforward evaluation of Eq. ( |3~T| ) using Eqs. ( j2.16j - |2.i5| ) leads to the result 

M 2 ^ 2 ' + UP- - PA 

(3,4|K|l,2) = 2(3,4|y|l,2> W== = f> , (3-5) 

V ™1 ^2 ^3 ^4 

where {Pi — Pf) = 5^ 2 \Pi± — Pf±)5(P i h — Pj~) and V is the standard expression for 
the sum of the 7r, and vector meson exchange potentials: 

(3,4|V|1,2) = (3,4\V{(f>) + V(tt) + V(V)\1,2). (3.6) 

The operator K is twice the usual two-nucleon potential times a factor which includes 
the light front phase space factor and a momentum-conserving delta function. Note that 
dk + /k + = dk 3 /E(k) for free nucleons where k + = E{k) + k 3 . 
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For the exchange of scalar and pseudoscalar mesons, only the term vxgo(P{ )v\ enters, 
and one finds 



n AWU ^ll 9\ ^(4)^(2) u(3)Tu(l) 

<3,4|V(*,7r)|l,2> = 4Af2(27r)3(g2 _^ 2) . (3-7) 
in which the momentum transfer q is given by 

q = k 3 - k x . (3.8) 

The notation is that u(i) is the Dirac-spinor for a free nucleon of quantum numbers i, and F is 
either of the form g s or i g^^T. The derivation of the contribution of vector meson exchange 
proceeds by including the meson exchange vxgo(Pi~)vx plus the meson instantaneous term 
t>3, and the result takes the familiar form: 

(3 4\V(V)\l 2) - 2 ^(4)7^(2)^(3)7^(1) (3Q) 
(3,4|V(V)|1,2) - 4M2(27r)3 ((?2 _ m?) • ( 3 - 9 ) 



The expressions (|3.7|) and (|3.9| ) represent the usual f27j-|3ll expressions for the chosen one- 
boson exchange potentials, if no form factor effects are included. The sum of the amplitudes 
arising from each of the individual one-boson exchange terms gives the invariant amplitude 
to second order in each of the coupling constants. The factors jtj? in Eqs. ( |3.7| ) and 



can be thought of as renormalizing the spinors so that uu = 1, and the factors J M: of 



k+ 

Eq. ( |3.5| ) serve to further change the normalization to u^u = 1 

These amplitudes are strong, so computing the nucleon-nucleon scattering amplitude 
and phase shifts requires including higher order terms. One may include a sum which gives 
unitarity by including all iterations of the scattering operator K through intermediate two- 
nucleon states: 

(3,4|M|1,2> = (3,4|K|1,2)+ £ / (3,4|iTl5,6) + + p ( P -\\ 7 ■ (5,6|A*|1,2). 

, , J p 5 p 6 r j — [p 5 -+- p G ) -t- le 



(3.10) 



The factor M f + l^jp-^ appears in both of the terms on the right-hand-side of Eq. (p.lOQ , 
so it is worthwhile to define a T-Matrix T using 

M* §{*,+) (p. _ p \ 

M ee 2T + V r ' (3 ' H) 

V ^1 ^2 ^3 ^4 

One realizes that Eq. (|3.10| ) is of the form of the Weinberg equation [Q (see Ref. ||) by 
expressing the plus-momentum variable in terms of a light-front momentum fraction a such 
that 

pt = aP?, (3.12) 
and using the relative and total momentum variables: 
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p_L = (1 - a)p 5± - ap 6± , 

Pi± = P5± + P6± ■ 



(3.13) 



Then, 

<3,4|T|1,2> = (3,4|V|1,2) + / £ (M^M) . <5,6|T|1,2), 



(3.14) 



where is the square of the total initial four-momentum, otherwise known as the invariant 

energy s and is the corresponding quantity for the intermediate state. Because the 

kernal V is itself invariant under Lorentz transformations in the three-direction and the 
integral involves p± and a the procedure of solving this equation gives T with the same 
invariance. 

Equation (1314]) can in turn be re-expressing as the Blankenbecler-Sugar (BbS) equation 



by using the variable transformation 



Ejp) + P 3 , . 

a = ^^T' (3 - 15) 



with E{p) = VP ' P + M' 2 . The result is: 

M 2 d 3 p 



(3,4|T|1,2) = (3,4|V|1,2> + / £ (3, 4|V|5, 6)— -- 2 ^—(5, 6|T|1, 2), (3.16) 

J E(p) p^ — p l + it 



which is the desired equation. Rotational invariance is manifestly obeyed. The three- 
dimensional propagator is exactly that of the BbS equation. There is, one difference between 
Eq. (|3.16|) and the standard BbS equation. Our one-boson exchange potentials depend on 
the square of the four momentum q 2 transferred when a meson is absorbed or emitted by a 
nucleon. Thus the energy difference between the initial and final on-shell nucleons is included 
and q° ^ 0. This non-zero value is a consequence of the invariance of D of Eq. ( |3.4| ) under 
Lorentz transformations in the three-direction. The usual derivation of the BbS equation 
from the Bethe-Salpeter equation specifies that q° = is used in the meson propagator. 
Including q° ^ instead of q° = increases the range of the potential relative to the usual 
treatment, and its consequences are explored below. One can convert Eq. (|3.16j) into the 



Lippman-Schwinger equation of non-relativistic scattering theory by removing the factor 



M/E(p) with a simple transformation [35 . 



B. Generation of a Realistic One-Boson Exchange Potential 

The present results are that one can use the light front technique to derive nucleon- 
nucleon potentials in the one-boson exchange (OBE) approximation and use these in an 
appropriate wave equation. Our purpose here is to show that the present procedure yields 



potentials essentially identical to the Bonn OBEP potentials p8| , p9| and these potentials 
lead to a good description of the NN data. 
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The Bonn one-boson exchange potentials employ six different mesons, namely, tt, T], u, p, a 
and the (isovector scalar) 5/a Q meson. The present formalism can account for the tt, 77,0; 
and a in an approximately chiral invariant manner. We wish to add in couplings ipr ■ Sip 
and ipT ■ p^'y^ip in a chiral invariant manner. Simply adding such terms to the Lagrangian 
of Eq. ( |2.1| ) would lead to a violation of the approximate symmetry of Eq. ( |2.3| ). However, 
one can redefine the operator U so that the symmetry remains. We replace the operator 
tjj'Uip' in the Lagrangian ( |2.1| ) by ^j'Uij)' : 

U = e^ T - p ^e^ 5 Ue^ T - 5 e^ T - p ^ . (3.17) 

Then the new Lagrangian is invariant under the transformation 

V/^e i75T 'V, U ->• e _il)5T-0 U e~* 7SJr, °. (3.18) 

In the present application we expand the exponential to first order in the meson fields. 

The final term we need to include is the tensor o pv q v part of the p-nucleon interaction. 
The presence of such a tensor interaction makes it difficult (or impossible) to write the 
equation for ip_ as ip_ = l/p + ■ ■ - ip + . This is a possible problem because the standard value 



of the ratio of the tensor to vector p-nucleon coupling f p /g p is 6.1, based upon Ref. |36 
Reproducing the observed values of E\ and P-wave wave phase shifts requires a large value 
f p /g p ; see Ref. [j37|. However our Lagrangian compensates for its lack of a p-N interaction 



with tensor coupling by generating such a term via vertex correction diagrams (which are 
the origin of the anomalous magnetic moment of the electron in QED). Such diagrams might 
not generate the phenomenologically required values of the coupling constants, but all that 
is needed here is that terms of the correct form be produced. This is because the standard 
procedure is to choose the values of the coupling constants so as to yield a good description 
of the NN scattering data. Thus we simply add in the necessary tensor terms. 

This brings us to the treatment of divergent terms in our procedure. The definition of 
any effective Lagrangian requires the specification of such a procedure. For the present, it is 
sufficient to say that we introduce form factors, F a (q 2 ) which reduce the strength of the a 
meson-nucleon coupling for large values of — q 2 . This is also the procedure of Refs. [55-50 . 
In principle, calculating the higher order terms using the correct Lagrangian can lead to 
consistent calculations of these form factors. We use a more phenomenological approach 
here. 

The net result is that the one-boson exchange treatment of the nucleon-nucleon poten- 
tial and the T-matrix resulting from its use in the BbS equation is essentially the same 



as the one-boson exchange procedure of Refs. p7L|30|]. The only difference is the keeping 
of the retardation effects — the square of the four-vector momentum transfer enters in our 
potentials. 



C. Specific One-Boson-Exchange Amplitudes 

The above formalism yields a one-boson-exchange potential (OBEP) which is a sum 
of one-particle-exchange amplitudes of certain bosons with given mass and coupling. Our 
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explicit expressions are presented here. As noted above, we use the six non-strange bosons 
with masses below 1 GeV/c 2 . Thus, 

V OB EP= £ Va BE (3-19) 

a=7v,rj,p,uj,ao,a 

with 7r and rj pseudoscalar (ps), a and ao/S scalar (s), and p and u> vector (v) particles. 

The OBE amplitudes (which are the contributions to V of our formalism) in the two- 
nucleon center-of-mass (cm.) frame are given by : 



(k'X 3 X i \V p ° BE \kX 1 X 2 ) 
9ps 



(2n)HM 2 
[k'X 3 X i \V s OBE \kX 1 X 2 ) 



u(k', X 3 )j 5 u(k, XiM-tf, X^u{-k, X 2 )[q 2 - m 2 ^]- 1 ; (3.20) 



(2tt) 3 4M 
[k'X^V^lkX.X,) 



a 2 

Js u(k',X 3 )u(k,X 1 )u(-k',X 4 )u(-k,X 2 )[q 2 -ml]' 1 ; (3.21) 



1 _ f v 

; {g v u(k', A 3 )7 At u(fe, Ai) + wh u ( k ', A 3 ) v(fc 3 - ki) v u(k, Ai)} 



(2tt) 3 4M 2 ^" v '-o/'f v 2M 
x{g v u(-k', X 4 )Yu(-k, A 2 ) + ^u{-k\ X A )a^i(h - k 2 ) v u{-k, A 2 )} 
x[q 2 -m 2 v ] 



2l—l 



1 _ f 

:{{9v + fv)u(k', A 3 )7 At u(fc, Ai) - ^Tju(k', X 3 ){k 3 + fei) A ,u(fe, Ai)} 



(2tt) 3 4M 2 v ' J/ v ' LJ 2M 

*{{9v + f v )u(—k', X^u(-k, A 2 ) - ^(-k', X,)(h + hru(-k, A 2 )} 
x^-m 2 ,]" 1 . (3.22) 
Our notation in the cm. frame is such that in-coming nucleon 1 carries helicity Ai and four- 



momentum k\ = k = (E, k) with E = \j M 2 + fc 2 , and in-coming nucleon 2 carries helicity 
A 2 and four-momentum k 2 = (E,—k); the out-going nucleons have A 3 , k 3 = k! = (E',k') 



with E' = \JM 2 + k' 2 , and A 4 , k 4 = (E f , -k') . The square of the four-momentum transfer 
between the two nucleons is q 2 = (k 3 — ki) 2 = {k' — k) 2 = (E' — E) 2 — (k' — k) 2 . The Gordon 
identity is used in the evaluation of the tensor coupling p9| . For the isospin- vector 
bosons Ti, a and p, the above amplitudes must be multiplied by T\ • r 2 . 

With an eye on the nuclear matter calculations to be conducted later in this paper, we 
note that in the factor f v /2M of the tensor coupling the nucleon mass M is used as a scaling 
mass to make the coupling constant /„ dimensionless. This scaling mass could be anything. 
Therefore, this M is not to be replaced by M* in the nuclear medium. 

In this subsection, we use Dirac spinors (in helicity representation) given by 

1 



u{k, Ax) = VE + M 2Al | fe | [Ax) (3.23) 

V E+M I 
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u(-k, a 2 ) = VeTm ( ) |A 2 ) ( 3 - 24 ) 

with 

I Ai> = x\i , |A 2 ) = X-a 2 , (3.25) 

where x denotes the conventional Pauli spinor. The normalization is u(k, X)u(k, A) = 2M. 
At each meson-nucleon vertex, a form factor is applied which has the analytic form 

/A 2 -m 2 \ na 

with m a the mass of the meson involved and A a the so-called cutoff mass; n a = 1 for 
pseudoscalar and scalar mesons and n a = 2 for vector mesons. Thus, the OBE amplitudes 
Eqs. (|3~20p-(|3T2l) are all multiplied by F 2 . 

D. Two-Nucleon Scattering 

In the two-nucleon cm. frame, the scattering amplitude T is the solution of the integral 
equation 

r M 2 1 

T(k',k) = V(k',k) + / d 3 P V(k',p) — - 2 T(p,k) (3.27) 

J hp k — p 2 + it 

where k, p and k' are the initial, intermediate and final relative momenta, respectively, 
of the two interacting nucleons; and E p = \/M 2 + p 2 . This is equation ( |3.16| ) with the 
spin indices suppressed for the purpose of simplicity. The corresponding equation for the 
i^-matrix (which we denote by R) is 



r M 2 1 

R(k',k) = V(k',k)+V / d 3 P V(k',p)—- 2 -R(p,k) (3.28) 

J £j v k — p l 



where V denotes the principal value integral. 

Using standard techniques |29| , |39| , the potential and the scattering equation are de- 



composed into partial waves. Numerical solutions are obtained by the matrix inversion 
method |29|j4l|] . For an uncoupled partial wave, phase shifts are then derived from the 
on-shell i^-matrix by 

tan 5\T lab ) = -~|fc|^r i* J (|fc|, |fc|) (3.29) 
with Tiab = Ik 2 jM and J the total angular momentum of the partial- wave state. For coupled 



partial waves and other technical details, see Ref. . 
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E. Results for the two-nucleon system 



Following established procedures [p8j , p9|l ) the coupling constants and cutoff masses of the 
six OBE amplitudes are varied within reasonable limits such as to reproduce the two-nucleon 
bound state (deuteron) and the two-nucleon scattering data below the inelastic threshold 
(about 300 MeV laboratory kinetic energy). 

In Table I, we show the meson parameters for our newly constructed Light-Front (LF) 
OBEP together with the predictions for the deuteron as well as low-energy neutron-proton 
(np) scattering. For comparison, we also give the parameters from an OBEP that was pre- 
viously constructed and applied in the Dirac-Brueckner approach to nuclear matter PRI30 



The latter uses the Thompson formalism |42| which is very similar to the BbS formalism — 
the propagator in Eq. ( 3.27 ) contains an extra factor of M/E p . Note that the Thompson 
OBEP uses n a = 1 also for vector meson form factors, which explains the differences in the 
vector meson cutoff masses between the two OBEP. 

Phase shifts for np scattering are shown in Fig. 1 for all partial waves with J < 2. Over- 
all, the reproduction of the NN data by our LF OBEP is quite satisfactory and certainly 
as good as by OBEP constructed within alternative relativistic frameworks. Based upon 
these results, we feel confident in applying this OBEP to the relativistic nuclear many-body 
problem. 



IV. NUCLEONIC TRUNCATION FOR THE MANY-BODY PROBLEM 

Now that the light front treatment of nucleon-nucleon scattering is in hand, we may 
proceed to the problem of computing the properties of infinite nuclear matter. We derive a 
light front Brueckner theory from first principles starting with the field-theoretic light-front 
Hamiltonian. 

The nuclear wave function for the ground state of infinite nuclear matter at rest is defined 
as 1^}, and we wish to solve the equation 

P-|*) = M A |¥), (4.1) 



in which P is the light-front Hamiltonian of Eqs. ( ^.Il[ )-( p.l8[ ). For a nuclear system at 
rest we must have also the result that 

P+\m) = M A \m). (4.2) 

It is necessary to discuss the light front Hamiltonian, and to find good approximate solutions 
of the above equations. 
We recall that 

P- = P (N) + J (4.3) 

in which Pq(N) is the kinetic contribution to the P~ operator, giving — ^ — for the minus- 
momentum of free fermions. The operator J is the sum of three terms of Eqs. ( |2.16|) -( p.l8|) : 



J = vi + v 2 + v 3 (4.4) 
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The operator v\ gives all of the single meson-nucleon vertex functions. The operator t>2 
accounts for instantaneous fermion exchanges: meson emission followed by instantaneous 
fermion propagation (propagator is ^+ ) followed by another meson emission. The operator 
v 3 accounts for the instantaneous propagation of vector mesons. 

We shall proceed towards an approximate solution of Eq. ( |4.1| ) , in two stages. We shall 
first consider the nucleons only part of the Hilbert space. This involves the assumption 
that using a nucleon-nucleon interaction K accounts for the meson-nucleon dynamics. This 
assumption is relaxed in Section VI, which displays the formalism necessary to construct the 
best possible potential and how to include meson degrees of freedom in the wave function. 
Our Hamiltonian (P~) contains no terms in which the vacuum can spontaneously emit 
particles. This simplifying feature causes the derivations to look very similar to those of 
non-relativistic theory, even though the treatment is relativistic. 



A. Introducing the two-nucleon force 

At present all of the interactions are expressed in terms of the meson-nucleon vertex 
functions and contact terms represented by the operator J. We shall follow the traditional 
path of using a two-nucleon potential and temporarily eliminate the meson degrees of free- 
dom. One way to accomplish this is to add and subtract the two-nucleon potential to the 
Hamiltonian and treat terms involving the difference between J and the two-nucleon po- 
tential as a perturbation. The use of light front dynamics mandates that we perform this 
operation on the Lagrangian because the construction of the Hamiltonian uses the field 
equations to identify the dynamical degrees of freedom, such as Therefore we need to 
study the effective Lagrangian: 

£ v = ip(ij-d-M)ip - —. (4.5) 

which removes the meson-nucleon interaction term J from the Lagrangian and replaces it 
with the density /C for the two-nucleon interaction K of the previous section. This means 
that 

K = i j cPx^dx-JC, (4.6) 

where K is given in Eq. (|3.5| ) and is twice the nucleon-nucleon potential times a kinematic 
factor. We recall that one must eliminate, using Eq. (|2.5| ), the components of interaction K, 
that connect ip_ to 

Given the new Lagrangian £y we may construct the corresponding P^ 1 operators using 
the canonical definition 

^—f^-mf*" (47) 

in which the degrees of freedom ip and ip are labelled by <p r . The term V involves only 
nucleon fields, not their derivatives, so the second term of the energy-momentum tensor do 
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not enter in computing T§ v . The element \iv = A — (g + = 2) is needed to construct the 
relevant P-minus operator and we find 

T+- = —2-0 (ij ■ d - M) tj) + K + 2iplid~ip + . (4.8) 

The origin of the factor | that multiplies K in the Lagrangian ( |4.5| ) is that /C enters here in 
T-jJ . It is worthwhile to define the P-minus operator obtained by using £y as P ~ with : 

p- = I J dx-d 2 x ± T+- = P -(iV) + K. (4.9) 

The complete P-minus operator is given by 

P- = P +H 1 , (4.10) 

with 

Pi EE J-K + p-( m ) (4.11) 

where P ~(m) accounts for the non- interacting mesonic contribution to P . The formal 
problem of choosing the best K by minimizing the effects of H\ = J — K is discussed in the 
Section VI. We shall assume here that the present OBEP is a reasonably satisfactory version 
of the best interaction, and we shall ignore the influence of the term H\ in calculations of 
the energy. 

The purely nucleonic part of the full wave function is defined as | $) , and is the solution 
of the light front Schroedinger equation 

P - | $) = (P -(iV) + K) I $) = M | $>. (4.12) 

The eigenvalue problem stated above is considerably simpler than the initial one, but does 
contain the full complications of the nuclear many-body problem. We shall next discuss the 
light front Hartree-Fock and Brueckner Hartree-Fock G5j approximations. 



B. Light Front Hartree Fock Approximation 



The equation ( |4.12 ) represents a difficult many-body problem. But the similarity between 



the light-front and equal-time results obtained for the nucleon-nucleon potential indicate that 
the same physical concepts are relevant, independent of the dynamical scheme. Therefore 
we use a scheme analogous to that of traditional Brueckner theory. The first step is to 
introduce a mean-field (MP) potential 



U 

2 



1 ( 

- J d 2 x±dx~ip(x) (U s (x) A- 7 • Py(x)) ip(x) 



- I d 2 x ± dx~i){x) lu s (x) + — U v {x) W(x), (4.13) 



a single-nucleon operator, with the second equation true for infinite nuclear matter in which 
the only non- vanishing component of Uy is U v = Uy = Uy = Uy . All quantities in the 
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integral of Eq. ( |4.13| ) are evaluated at the same value of x + , chosen to be 0. The operator U 
is to be determined ultimately by the light-front G-matrix defined below. The idea is that 
U can be chosen so as to provide a good representation of the effects of the two-nucleon 
interaction K. The mean field Lagrangian Cmf is defined by removing the effects of K and 
replacing these by the effects of U. Therefore we may specify 

C MF = i; (z 7/1 3" - M)ip - $(x) (u s (x) + ^-Uyixij ^(x). (4.14) 

This Lagrangian leads to the nucleon field equation: 

(id- - Uy) = (a ± ■ P± + P(M + U s ))il>- 

id+il>_ = (<x ± ■ p ± + (3(M + U s ))1> + , (4.15) 



in which we have made the mean-field version of the transformation ( |2.5| ) with 

d+g v A M F = U+{x). (4.16) 

The light-front Hamiltonian density Ty~~ can now be obtained from Eq. ( |4.8[ ) using the 
field equation ( |4.15| ) as 

T+- = 2^ + id-tP + - 2^ (^-U v + U s ]ip + K. (4.17) 

It is also worthwhile to obtain the plus-momentum density Ty" + which is 

T++ = 24> + id + i) + . (4.18) 

The purpose of introducing the mean field approximation is that the eigenvalues and 
eigenvectors of light-front mean field Hamiltonian P^ F are easy to obtain and can be chosen 
so as best approximate the effects of the two-nucleon interaction. The mean field light front 
Hamiltonian density T^ IF is obtained from Eq. ( (4.7|) using Eq. ( f4.14|) (or by setting /C to 
zero in Eq.( |4.17| )) as 

T+ F = 2^ + icTV>+ - 2^ (^—Uy + U s ^j V, (4.19) 
and its volume integral is Pmf'- 

Pmf = \J d 2 x ± dx-T+ F , (4.20) 
a single-nucleon operator. Setting K to zero in Eq. ( |4.17|) and using Eq. (|4.9|) shows that 

Pmf = P (N). (4.21) 
The ground state eigenvector of this operator is a Slater determinant denoted as | (j>): 

Pmf I <P) = Pq(N) \ <P) = m \ 0). (4.22) 
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We shall use both the Hartree-Fock and Bruckner Hartree Fock approximations to obtain 
expressions for U. For now we pursue the question: Given a U, how do we proceed? The 
first step is to expand the field operator ip in terms of the eigenfunctions of id~ in the light 
front Dirac equation ( |4.15|) The nucleon field operator is constructed as follows: 



= r d 2 k ± dk+6(k+) y, k . 
V ; J (2tt) 3 / 2 v / 2F : a 

where k-x = \ {k~x + + k + x~) — • a?j_. We keep only the nucleon part of ip(x) as the anti- 
nucleon degrees of freedom are not needed here. These spinors u(k, A) are the eigenfunctions 
of Eq. ( |4.15| ), with normalization u(k, X)'~f + u(k, A) = 2k + . For the present treatment of the 
translationally invariant infinite nuclear matter system, the mean-field potentials Us and 
Uy are independent of the spatial position x±,x~. The eigenvalues of Eq. ( [4.15[ ) are given 

by 1,| 

fcj + (M + Usf 
k —U v -\ — , (4.24) 

in which Us and Uy depend upon k± and k + . 

The next step is to better define the Slater determinant | <fi). The occupied states are to 
fill up a Fermi sea, which is usually defined in terms of a Fermi momentum, kp, that is the 
magnitude of a three vector. This three vector is defined as: 



k + — y (M + Us) 2 + k- k + k , (4.25) 
which implicitly defines k 3 . Using Eq. ([4.25|) allows one to maintain the equivalence between 



energies computed in the light front and equal time formulations of scalar field theories |47 
and to restore manifest rotational invariance in light-front QED |48| . 



The computation of the energy and plus momentum distribution proceeds from taking 
the appropriate expectation values of the energy momentum tensor Ty U . 

P^ = \J d 2 x ± dx-(<P | T+" | 0). (4.26) 

We are concerned with the light front energy P~ and momentum P + . The relevant compo- 
nents of Ty U are presented in Eqs. ( 4.17 ) and ( ^.18| ). Taking the nuclear matter expectation 



value of Ty and Ty + and performing the spatial integral of Eq. ( |4.26| ) leads to the result 
P v- 4 ' rfW( fci + ( f + +C/s)2 - 2 l£ gftj)^^} +{4 ,\K\4,) (4.27) 



!i (2i) 3 Jf [ k+ 2 Y v^Ff -M 

where fl is the volume of the system Q = | / d 2 x±dx~ . The subscript F denotes that 
\ k \< kp with k 3 defined by the relation fl4.25|) . The integral involving u(k, X)Usu(k, A) 
may also be expressed as: 
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Us 



(27T) 



d k i dk~* 



1 ^ u(k, A) A) 



(4.29) 



Equations (|4.27|) and ( 4.28|) along with the expression for k + , ( 4.25 ) allow an evaluation 
of P~ and P + . This shall be done by obtaining the mass M of the A-nucleon system 
as M = \{Py + -Py ) and minimizing Mq per nucleon. For a nuclear system at rest, the 
eigenvalues of the plus and minus momentum operators must be the same, so it should be 
legitimate to use the average value. In practise, approximations are made. However, the 
procedure of minimizing the average \{Py + Py ) is the same as minimizing the expectation 
value of Py subject to the constraint that the expectation value of Py is equal to the value 
of Py [0 . Summing equations ( |4.27|) and ( |4.28| ) and dividing by a factor of 2 leads to 



Mo 



4 1 



(2tt) 3 2 



d 2 k , dtf 



'k 2 ± + (M + U s 

k+ 



+ 



K 

y 



(4.30) 



Then replace the integration over k + by one over k 3 using Eq. ( |4.25|) and the definition 

El = y/k-k + (M + U s y, (4.31) 
so that Eq. Q4.30D takes the form: 
M 



with 



(2?r) 



U, 



J d 3 k6{k F - k) El - (<p I U s I <p) + 



K 
~2 



(4.32) 



(27T) 



d 3 ke(ki 



2Ei 



2Et 



(27T) 



2E* k 



(4.33) 



One obtains a formalism that looks more conventional by using a discrete representation 
of the single nucleon states. We define a set of spinors | a), with a representing the quantum 
numbers k and A such that 



U' 



a) 



■u(k, A) 1 



2E* k 



(a 



{a | 7 
1 
2 







d 2 xj_dx (a \ x)(x \ a). 



The difference between (x 



(4.34) 
cx)eTi of the 



a) and a usual equal time (ET) spinor, (x 
same quantum numbers, energy e(k) = El + {7y, and normalization can be determined by 
considering the phase factor, using Eq. (|4.24|) : 



k ■ x 



x 1 



+ 



k+ 



x 



k 



X 



Uv + 



kl + (M + U s ) 2 \ (t + z) , ,,(t-z) 



k+ 



• x 



e(k)t — k ■ f — i 



UW 



(4.35) 
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The last factor is the consequence of using the barred form of the vector potential according 



to Eq. ( 4.16 ). The only difference between the light front spinors and those of the equal 
time form is due to this phase factor. The consequences of this phase factor for computa- 
tions of the light front momentum density are that nucleons carry only 65% of the nuclear 
plus momentum in the mean field calculation A similar result is to be found below. 

Furthermore this phase factor has the desirable feature of suppressing the number of nuclear 
anti- nucleons |4"9"|] . 

For now, we consider the effects of using light front spinors on the calculation of the 
energy. We express the energy M of Eq. ( ^4.32| ) in terms of the spinors | a) of Eq. ( |4.34j ) as: 



in which | af3) a =| a(3)- 



5> 

a<F 



a,/3<F 



(5a). In this notation 

E a = e a - (a \ 7°[/£ | a). 



(4.36) 



(4.37) 



so that we obtain 



a<F 



5> 

a<F 



U S + 1°U° V 



a ) + \ £ ( a P 

a,(3<F 



V I a(3) a . 



(4.38) 



This result gives the value of the nuclear energy in terms of the eigenvalues of the light front 
Dirac equation and in terms of light front spinors. The term K/2 has been replaced by V 
according to Eq. (|3.5|) and using the normalization of Eq. (^4.34|) . 



It is necessary however, to consider the effects of the phase factor 



of Eq. (EL35D 



which accounts for the difference between light front and equal time spinors. This factor 
and its complex conjugate multiply to unity in the calculation of the matrix elements (a \ 
Us + j°Uy | a) and in the matrix element (a/3 \ V | a,/3). We need to consider also the 
matrix element {a(3 \ V \ (3a) for which the effects of the phase factor do not automatically 
cancel. In principle, Uy is a function of the momentum denoted by the quantum numbers 
a, p. In practise, this dependence is weak and can be ignored in calculations of the energy, 
provided one evaluates the potential at a reasonably chosen average value. Furthermore, 
in evaluating the matrix element {a/3 \ V \ (3a) both of the states (3a are below the Fermi 
sea and have a momentum separated by an amount small compared to the scale of the 
momentum dependence. Thus the phase factor does not enter in present calculations of the 
energy (but does in the evaluation of the plus-momentum distribution). Our calculations of 
the energy in the light front and equal time formulations yield the same results. However 
calculations of the plus-momentum distributions can only be done using the light front 
formalism. 

Let us determine the mean field, and the corresponding value of M in terms of V. 
The light front Hartree-Fock (HF) approximation is defined by taking the mean field to be 
calculated from the average potential according to 



(a 



(TT i ^,0tt0\HF 

{Us + 7 U v ) 



a) = U(a) 



{a? | V | a(3) a . 



(4.39) 



I3<F 
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Summing over the occupied orbitals gives 

E^(«) = 2(0| V |0). (4.40) 

a<F 

In this case, the use of the above equation in the expression ( |4.38| ) for M , leads to the HF 
approximation for the nuclear energy: 

M HF = E ^ - \ E I V I afla, (4-41) 

a<F A a,/3<F 

which has the same form as the expression for the energy in the usual equal time Hartree 
Fock expression. 

We shall also need to obtain m , the eigenvalue of Pj^p. This expression will be used in 
obtaining the Bruckner Hartree- Fock approximation. According to Eqs. (|4.17|) and ( [4.22[ ), 
the difference between Mq and mo is the expectation value of the potential V. Thus 

™o = E ^ - E (« I U S + 1°U° V | a) (4.42) 

a<F a<F 



C. Light Front Brueckner Hartree-Fock Approximation 

The interaction K between two nucleons is strong and the scattering amplitude is ob- 
tained, as discussed in Section 3, by solving the Weinberg equation — the light front version 
of the Lippmann-Schwinger equation. Thus we need to go beyond the Hartree-Fock approx- 
imation. This shall be accomplished by treating the interaction between two nucleons to all 
orders in K. 

The idea is that we wish to find the Slater determinant | 0), recall Eq. ( [4.22j ), that leads 



to the best approximation for the energy M$ of the full nucleonic wave function | $), recall 
Eq. O. 



Both of the states | 0) and | $) are eigenstates of a P-minus operator, and both are 
eigenstates of the operator P + (A). We shall use standard techniques to derive a perturbation 
theory in K to obtain an expression for the state | $) in terms of | 0). Thus we write 

| $} = | 0) + A | $} (4.43) 

with 

A = 1-|0)(0|. (4.44) 
Then use Eq. ( |4.43|) in Eq. ( |4.12| ) and multiply the result on the left by A to obtain 



A 1 *> = M -AK(N) + K)A AK I < 4 - 45 > 



so that 
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M -A{Pq(N)+K)A 



AK 



(4.46) 



We can obtain a useful expression for M by acting with the operator (0 | (P (N) + K^j on 
the left of Eq. ( 4.46|) and using the result {(f) | $) = 1, which follows from Eq. fl4.46j ). Then 



we find 



M 



Pn+K 



+ 



KA 



M -A(P (N) + K)A 



AK 



Using Eq. ( ^22|) in Eq. (pL47p leads to 



M - m = (0 | K 
which can be re-stated as 



+ 



KA 



M - m 



where 



X = K + KA 



M -A(Pq(N) + K)A 



X 



AX. 



AK 



M Q - APq(N)A 



(4.47) 



(4.48) 



(4.49) 



(4.50) 



The operator X is a many-body operator acting on all nucleons via the iterations of the 
two-nucleon interaction K. We shall make the independent pair approximation of including 
only pair-wise interactions. Thus we approximate 



X 



(4.51) 



where is a two-nucleon operator which is a solution of the integral equation 

A 



Ka + Ki 



P,7-AP -(iV)A 



r«(P<7). 



(4.52) 



The notation i, j refers to a pair of particles. The relevant matrix element is expressed using 
the eigenstates of Eq. ( |4.15| ) as 



(3,4|r(P 1 ; 2 )|l,2) = (3,4|ir|l,2)+ ]T |(3,4|tf|5,6) 



2M* 



d 2 p 5± dp£ Q 



Pt Pe p i,2 - (P5 + P% ) + * e 



in which we define 



M* = M + C/g. 



(5,6|r|l,2), 
(4.53) 



(4.54) 



The operator Q, to be specified below, is the two-body version of A and projects the momenta 
P5 and p6 above the Fermi sea. The factor — f + ^ P ^~^ appears in each of the terms of 

rC-^ ^2 ^3 ^4 

Eq. ( [4.53[ ), so it is worthwhile to define a Bruckner G-Matrix G using 
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M* 2 5^ +) (P- - PA 
r = 2G -. 1 1 f) . (4.55) 

'h + h + h + h + 

A/^ rL<2 ^3 ft- 4 



To follow the steps of Sect. Ill in converting Eq. ( [4.53 ) into one of a more familiar form, 



in which rotational invariance is manifest, one needs to know the values of 

P^=Pl+P%, ( 4 -56) 

which for the case of relevance here in computing the nuclear expectation value in the inde- 
pendent pair approximation, is the same as p§ + p$ ■ The single-particle minus- momentum 
eigenvalues are given according to Eq. ( [4.24]) as 



Our approximation is that Uy is independent of orbital i. Thus this potential cancels out 
in computing the difference P[~ — (p$ + p$) and the energy denominator is as in the free 
space considerations of Sect. Ill, except that the mass of the nucleon is replaced by M + Us- 
Thus the previous derivation of an equivalent three-dimensional integral equation that is 
manifestly covariant and rotationally invariant proceeds as before. 

One expresses the plus-momentum variable in terms of a light-front momentum fraction 
a of Eqs. (|3.15|) and ( |3.13|) so that one obtains 



r 2M* 2 d 2 r)\dotG 
(3,4|Gil,2> = (MV|l,2}+/ £(3,4|V|5,6>— /^"j 2 (5, 6|G|1, 2), (4.58) 

where Pf is square of the total initial four-momentum, computed using M + Us for the 
nucleon mass. 

Equation ( |4.5^ ) can in turn be re-expressed as a medium-modified Blankenbecler-Sugar 



(BbS) equation |33| by using the medium-modified version of the variable transformation 



E* + v 3 

a = -JLL*-, (4.59) 
2E* v ' 



with E* given in Eq. (|4.31|) . The result is: 

r M* 2 d 3 n O 

(3,4|G|1,2) = (3,4|\/|1,2)+ / £(3,4^5,6) 2 P \ (5, 6|G|1, 2), (4.60) 

A 5 ,A 6 Pi P 

which is the desired equation (with Dirac spinors normalized as in Sect. III). 

The Brueckner light front Hartree-Fock (BHF) approximation is defined by taking the 
mean field to be calculated from the average G-matrix according to 

U(a) = (a | {U S + i Q U° v ) BHF \ a) = £ (a/9 | G \ a(3) a . (4.61) 

I3<F 
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where Dirac spinors are normalized as in Eq. ( |4.34| ). The sum over occupied orbitals gives 



G 



(4.62) 



We use this BHF mean-field to determine the value of m via Eq. ( [4.42 ). Then the use of 
Eq. ( }4.51|) in Eq. ( |4.49|) determines the value of M as the eigenvalue of Py . But M is 
also the eigenvalue (or in this case the expectation value) of Py. The minimization of Py 
subject to the constraint that the expectation value of Py~ is the value of Py leads to the 
BHF version of Eqs. ((OT|) and flO§: 



P + 



(2k) 



d 2 k<dk A 



kl + (M + U s 



k+ 



1 y, u(k, X) LT u(k,X) 

2 Y v 7 ^ 5 v 7 ^ 



n (2vr) 3 7f 



d 2 k l( ik + k + . 



(4.63) 
(4.64) 



Note that the quantity k + is defined in Eq. ( |4.25|) . 

Taking the average of equations ( |4.63| ) and ( [4.64| ), and using the basis of Eq. (|4.34|) leads 
to our result for the BHF version of the nuclear mass: 



a<F 



E ( a P I G | a/3) , 



(4.65) 



a,f3<F 



This is equivalent to the usual expression of Ref. [p8lj3~0|l , cf. Eqs. (|5.10|) and ( |5.11|) below. 



V. LIGHT-FRONT BRUECKNER THEORY OF NUCLEAR MATTER 



A. Summary 

The formalism of the previous section can be summarized using the notation of 



Ref. [28,30]. In that work, single- nucleon motion in nuclear matter is described by the 
Dirac equation 



— M — U)u*(k, s) = 



or in Hamiltonian form 



[a ■ k + f3M + pU)u*(k, s) = e k u*(k, s) 



(5.1) 



(5.2) 



with 



U = U S + 1 °U° V] (5.3) 
where we use the usual notation p8| , (3 = 7 and a = 7°7. The solution of Eq. (|5.1|) is 



Ei + M* ( 1 \ 



IE, 



fc \ Et+M* 



(5.4) 
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with M* defined in Eq. ( |4.54j ). E%, given by Eq. ( }4.31|) , and Xs a Pauli spinor. The normal- 
ization is u*(k, s)u*(k, s) = and u*\k, s)u*(k, s) = 1, as in Eq. ( |4.34| ). Notice that the 



in-medium Dirac spinor Eq. ( p.4| ) is obtained from the free Dirac spinor by simply replacing 
M by M*. The difference between the equal time spinors used in Refs. [|28|j30|1 and light front 
spinors is the phase factor discussed above in Sect. IV. C. As noted, this factor cancels out of 
the matrix elements when using the approximation that the vector potential is independent 
of momentum. 

The nuclear matter G-matrix is the solution of the integral equation, ( |4.60| ) without the 
explicit spin indices, 

G(k', k- \p\,h F ) = v*(k', k) + / ^-v*{k',p) M * 2 y \p\M , 

(5.5) 

where P is the total three-momentum of the two nucleons in the nuclear matter rest frame; 
and k, p and k' are the initial, intermediate and final relative momenta, respectively, of 
the two nucleons interacting in nuclear matter. Note that the third components of the 



various momenta are obtained by using the appropriate versions of Eq. ( |4.25| ). kp denotes 
the magnitude of the Fermi momentum corresponding to the nuclear matter density under 
consideration. The Pauli operator Q, for which we use the angle average Q ||41|| , projects onto 
unoccupied states. In the derivation, we have used the angle averages (|P±fc) 2 « ^P 2 + k 2 

and ( \P ± p f ss \P 2 + p 2 . The latter implies E\ p+p « ^ M* 2 + \P 2 +p 2 . 

The essential difference with standard Brueckner theory is the use of the potential V* in 
Eq. (|5.5|) . As indicated by the asterix, the OBE potential of Eqs. ( |3.20|) -( |3~2~2"1) is evaluated 
by using the in-medium spinors Eq. (5T4") instead of the free ones. 



It is necessary to discuss a technical point concerning the retardation effects in the form 
factors Eq. ( |3.26| ) that enter in evaluating V*. Consider the matrix element of V* for which 
our formalism says to use 

q° = E* a - E; = (E* a + U v (a)) - (E* p + U v (a)). (5.6) 

We have consistently ignored the state dependence of Uy. This is a good approximation for 
energy differences governed by the Fermi energy. However, the relevant energy scale in the 
form factors are the parameters A of Table 1, which are on the order of a GeV or more. 
Thus for these terms it would be better to use 

q° = {E* a + Uy(a)) - (E; + U V (P)). (5.7) 

The quantity E^ + Uy(a) for occupied orbitals a is close to the nucleon mass. For high 
energy orbitals (3, U{(3) is small. Thus the use of Eq. ( |5.7|) is numerically very similar to 
ignoring the effects of the medium modifications in evaluating the retardation effects in the 
form factors, which is what we do. 

The single-particle potential for nucleons in nuclear matter is 

M* 

U(a) = (a\U\a) = (a\U s + l°U v \a) = —U s + U v ; (5.8) 



24 



which is calculated from the G-matrix by 



U(a) = Re J2(aP\G\aP) a (5.9) 

where a denotes a state below or above the Fermi surface (continuous choice). 

Using the notation of the previous section, we define 'the energy per nucleon in nuclear 
matter' by 

which can also be written as 

-a = T E H7 ■ K + M\a) + t^t E (*P\G\<*P)a - M . (5.11) 

a<F a,f3<F 

Note that this equation depends on the Fermi momentum kp and, thus, on the density of 
nuclear matter. 

It is useful to have the following summary of formulae concerning the single particle 
energy e a : 

e a = (a\j ■ k a + M + U\a) (5.12) 
= (a|7 • k a + M\a) + U(a) (5.13) 

= ™£*+£u, + U} (5,4) 

a a 



K + K, (5-15) 



with E* = JM* 2 + k 2 n and M* = M + U s . 



The calculation of the nuclear matter G matrix involves a self-consistency, since the 
solution of Eq. ( |5.5| ) for G requires knowledge of M* which, in turn, is determined from G 
via Eqs. (|5.8| ) and ( |5.9|) . In practise, one starts out with an educated guess for M*, solves 
Eq. (|5.5| ) for G and uses this G to calculate a new M* from Eqs. ( |5.8| ) and (|5.9| ). The 
procedure is then repeated starting with the new M*. This is reiterated until the calculated 
M* reproduces accurately the starting M*. 



B. Results 

The formalism of the previous section is used to calculate the energy per nucleon in 
nuclear matter as a function of density, Eq. ( |5.11| ). Our result is plotted in Fig. 2 by the 
solid line. The curve saturates at £ /A = —14.71 MeV and kp = 1.37 fm _1 , and predicts an 
incompressibility of K — 180 MeV at the minimum. These predictions agree well with the 
empirical values S/A = -16 ± 1 MeV, kp = 1.35 ± 0.05 far 1 , and K = 210 ± 30 MeV |0 



To get a better idea of the quality of our predictions, it is useful to compare with 
the results from alternative relativistic approaches. Brockmann and Machleidt |3(J predict 



S/A = —13.6 MeV, kp = 1.37 fm , and K = 250 MeV at saturation, using the equal-time 
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formalism and their 'Potential B'. The greatest difference occurs for the incompressibility 
which is predicted smaller by the LF Brueckner theory implying a softer equation of state. 
This can be partially attributed to the medium effect that comes from the meson propagators 
in the LF approach and that is absent in the equal time (ET) approach. Recall that in the 
LF formalism the momentum transfer between two nucleons exchanging a meson is 

q=(q ,q) = (E'-E,k / -k), (5.16) 

where E and E' are nucleon on-mass-shell energies [more explanations can be found below 
Eq ( p.22j )1, implying the meson propagators 



(5.17) 



q 2 - m 2 a (E' - E) 2 - (k' - k) 2 - m 2 a ' 
while in the ET formalism no energy is transfered, thus, 

q = (0,q) = (0,k'-k), (5.18) 



and the propagators are 



—q 2 — ml —{k' — k) 2 — ml 
In nuclear matter, the free-space LF meson propagators Eq. ( [5.17 ) are replaced by 



(E'* - E*) 2 - (k 1 -k) 2 -m 



(5.19) 



(5.20) 



while the ET propagators undergo no changes. The medium effect on the LF meson propa- 
gators enhances them off-shell which leads to more binding energy. This is demonstrated in 
Fig. 2 where the difference between the dotted and solid curve is generated by the medium 
effect on the meson propagators. 

There is another difference that arises from a technical issue in the solution of the 
transcendental equation for the G-matrix. We obtain new values of the mean fields 
M*(k F ) = M N + U S = 718 MeV and U v = 165 MeV. The mean field potentials obtained 
here from the G-matrix are considerably smaller than those of mean field theory in which 
the potential is used. We discuss the implications for nuclear deep inelastic scattering in 
Sec. VII. 

The most important medium effect in relativistic approaches to nuclear matter comes 
from the use of in- medium Dirac spinors representing the nucleons in nuclear matter ('Dirac 
effect'). This effect (and the medium effect on meson propagators) is absent in the con- 
ventional (nonrelativistic) Brueckner calculation which yields the dashed curve in Fig. 2. 
Characteristic for all predictions using conventional Brueckner theory is that the saturation 
density is predicted too high and, thus, they all fail to explain nuclear saturation correctly. 

The effect that is generated by the in-medium Dirac spinors is strongly density dependent 
(due to the density dependence of M*) shifting the saturation curve towards lower densities 
such that nuclear saturation is predicted at the correct density (solid and dotted curves in 
Fig. 2). 
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The effect from the in-medium Dirac spinors is, of course, largest for a and u exchange 
for which the LF and ET formalisms predict essentially the same. However, ir and p do also 
contribute to the medium effect and, here, we have differences between LF and ET. The 
general underlying reason for this difference is that for derivative coupling, implying a mo- 
mentum dependence of the meson-nucleon vertex, the difference in the momentum transfer 
(= meson momentum) between LF and ET [Eqs. (|5.16 ) and ( |5.18|) , above] creates a differ- 
ence in the vertices (and OBE amplitudes) — besides the one in the meson propagators. The 
p includes the tensor coupling i-^cr^ u q v and the LF OBE amplitude is given in Eq. ( |3.22| ). 
However, in the ET formalism, the p amplitude consists of Eq. ( |3.22| ) plus additional terms 
that contribute only off-shell, as discussed |39|]. In nuclear matter, the medium effect gen- 
erated by these off-shell terms essentially cancels the medium effect that comes from the 
main part of the amplitude, Eq. ( |3.22| ). Therefore, ET p exchange produces a much weaker 
medium effect than LF. 

Concerning the pion, the pseudovector (pv) coupling (or gradient coupling) i ^ N 757 M g M 
has been generally used, in the ET formalism ||30| . The resulting one-pion-exchange (OPE) 
amplitude can be cast into a form that consists of the amplitude Eq. ( |3.20| ) plus off-shell 
terms []51) , [52]] . In the LF formalism, no matter if the pv or ps coupling is used, the OPE 
amplitude always comes out to be Eq. ( p.2G|) , and there are no additional terms |p2|| . In 
nuclear matter, the medium effect from the off-shell terms of the ET formalism damp the 
medium effect from the main OPE amplitude, Eq. (|3.20|) , similar to what happens with the 
p. Therefore, LF ir exchange produces a stronger (more repulsive) medium effect than ET. 

In summary, the Dirac effect comes out more repulsive in the LF formalism as compared 
to ET due to off-shell differences in the 7r and p exchange amplitudes. On the other hand, 
the LF formalism generates an attractive meson propagator effect that is absent in ET. As 
it turns out, these two effects cancel to a large extent, leading to a non-trivial similarity 
between the LF and ET results. 



VI. FULL WAVE FUNCTION AND MESON DEGREES OF FREEDOM 



The nucleonic wave function | $) has been determined in Eq. ( f4.46| ). This gives us the 



purely nucleonic part of the Fock space, in which the effects of the mesons have been replaced 
by the two-nucleon interaction K. However, the full wave function is | of Eq. ( f4.1| ). We 
need to assess whether | <3>) is a good approximation to | and we also need to determine 
the mesonic plus-momentum distributions. 

We recall the relation between the full P~ operator and the one of Eq. ( |4.12| ) (P ~ = 
Pq{N) + K) corresponding to the nucleonic wave function | $): 

P~ = P +J-K + P (m). (6.1) 

Using this in Eqs. ( [4.1|) and ( |4.12j ) allows us to obtain 



* >H *> V< - aVa/ * (J ~ 70 (6 ' 2) 



where A$ = 1— | <&}($ 
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An expression for the nuclear mass M A can be obtained by multiplying Eq. ( |6.2| ) by 
($ | P~ using ($ | = 1 to obtain the result: 

M A = ($ | P- | $) + ($ | (J - K)A $ __J-^- x -A $ (J - K) | $). (6.3) 

For the purely nucleonic wave function | $) we have 

($ | p~ | $) = M + ($ | J-K | (6.4) 

so that 

M A = M + ($ | J - K | <£> + ($ | ( J - A)A $ — 1 — A„(J - A") | $). (6.5) 

The difference between M A and M is the expectation value of the operator O: 

Q S J-KHJ- *)^ Ma _ Wt \ At(J _ K)A MJ - K), (6.6) 

and a bit of algebra shows that O satisfies the integral equation 

O — J — K+(J- K)A< 5 >G (M A )A®O (6.7) 

where 

G (M A ) = A*- 1 A«. (6.8) 

Mj4 — i\.$-Tg A<j> 



The lowest relevant order of Eq. ( |6.7| ) is given by 

O^J-K + (J- K)A^G (M A )(J - K). 

The one-boson exchange interaction K is given by Eq. (|3.1| ) (which now includes also 
the effects of Sect. III.B) and we may determine if the expectation value (<3> | O | $) is 
reasonably small. If this is true, the quantity Mq would be a good approximation to the 
true eigenvalue of the P~ operator M A . In the one-boson exchange approximation 

($ | O | $) « ($ | v 3 - K | $) + ($ | J G (M A )J | $} (6.9) 
= ($ | w 3 - K | $) + ($ | uiG (M A )ui | $). (6.10) 

The use of Eq. Q3.1| ) yields 

($ | O | $) « ($ | «! (G (M a ) - ^(^7)) «i I (6.H) 

where the term Pj~ is specified in Eqs. ( |4.53| )-( f4.57| ). But within the independent pair 
approximation (in which one includes only the energy (minus-momentum) differences for a 
chosen pair of nucleons) 

G (M A )=g (Pr), (6.12) 

and the expectation value of O vanishes. 

Thus within our approximations, it is consistent to say that the exact nuclear mass M A 
is well approximated by M . This means that we have shown that it is acceptable to remove 
the explicit mesons for calculations of the nuclear mass. One could evaluate the presumably 
small corrections by going beyond the independent pair approximation. The simplicity of 
the derivation of this result is made possible by the dynamical simplicity of the vacuum, 
which is one of the defining features of light front field theory. 
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A. Momentum distributions 



We can now compute the + components of momentum. Look at T ++ as given by 
Eq. ( |2.19|) The plus momentum carried by the scalar meson is given by 

p+(«/>) = J d 2 k ± dk + k + a j {k)a{k), (6.13) 

while that of the pion is given by 

P+(tt) = J d 2 k L dk + k + a)(k) ■ a{k), (6.14) 

and that of the vector meson is given by 

p+( w )=^ [ d 2 k ± dk + k + a\k,u)a{k,uj). (6.15) 

w=l,3 

We shall handle the scalar term first. The evaluation of Q6.13Q using (|6.2| ) leads to 

P + (0) = J d 2 k ± dk + k+(<S> | j(fc) — — -1^— j(fc) | $) 

« / d 2 k ± dk + k + ($ | i(fc) (MA _^ p _ A)2 i(fc) | *>, (6.16) 

in which the approximation is motivated by the near equality of Ma and Mq and the universal 
expectation that the impulse approximation evaluation of the meson exchange potential is 
valid. The term j(k) is defined via the contribution of the scalar mesons to V\. 

Vl ((f)) = J d 2 k L dk + [j(k)a(k) + H.C], (6.17) 

in which j(k) can be obtained from Eq. fl2.16|) and is a nucleonic operator that depends on 
k± and k + . 

The operator to be evaluated in the above equation has one and two body pieces. The 
one-body terms are related to a shift in the self energy of the nucleon caused by the medium. 
In infinite nuclear matter the ratio of pairs to single nucleons is infinite so that the number 
density is well approximated by the two nucleon terms of Eq. ( |6.16| ). The evaluation is 
simplified by the use of Eq. ( |6.12|) and noting that the relevant matrix element is the same 
as occuring in the one-boson exchange operator K except that the denominator is squared. 
Thus the momentum density n^(fe), defined by 

P+(0) = j d 2 k ± dk + k + n<t) {k), (6.18) 

is given as a derivative of the scalar-meson exchange contribution to the nucleon-nucleon 
potential: 

n,(k) « -2($ | ^ (fc)) | $) (6.19) 
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with 



avW7( fe )) 

<9P,7 



3i(k) 



(M A - AP-A) 2 



3j(k) 



(6.20) 



in which the notation i,j specifies that only two-nucleon contributions are included. Note 
that | $ > is the correlated ground state. Note that the expectation value is taken using 
the single particle basis specified by Eq. (|4.34f) . We recall Eqs. ( |3.1|) and (|3.2j ), and use 



k + (P^-P )=q 2 -ml. 
Note that the momentum of the exchanged meson is k, and 

k 2 + m 2 



k 



q , k_L = q_L,k~ 



k+ 



where q is the nucleon momentum transfer. Thus one may obtain the result that 

d 1 k + 



2\2 ' 



(6.21) 



(6.22) 



(6.23) 



This means that evaluating the plus-momentum distribution for scalar mesons is the same as 
evaluating the expression for the scalar meson contribution to the nuclear potential energy, 
except that the potential is multiplied by the factor —k + /(q 2 — m 2 ,). The net result is that 



n,(fc) = £ (a/3 | n{ p V^(k) 



-k- 



a,/3<F 



(q 2 - m 2 ) 



(6.24) 



where £l a ,i3 is the Moeller scattering operator for the two-nucleon state a(3. One finds a 
similiar expression for the pionic density with 



7T (k)= £ (a/3 | fl^K(fc) [ _ fe+) 2 ^ a/3 I a/9),. 



(6.25) 



a,f3<F 



The evaluation of the vector meson density n^(k) requires more steps because the meson- 
exchange potential has a contribution from the instantaneous meson exchange term i> 3 . This 
instantaneous term does not lead to a meson "in the air" and therefore does not contribute 
to n U) (k). One finds that 



n 



,{k) = £ (a{3 | nlpV u (k) 



-k + ) 



a,/3<F 



(q 2 - ml 



-Vt a p I aft) a, 



where 



V u (k) = g 2 J(2nf F * {q2) 



(q 2 — ml) 



k+< 



(6.26) 



(6.27) 



It is also necessary to discuss the nucleonic plus-momentum distribution. This is deter- 
mined in Ref. 0. Here the nucleon-nucleon correlations cause the momentum density to 
have contributions from above the Fermi sea. We have 
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Ptr 



N 



A p B (2n) 



J d 2 k ± dk + k + N{k ± ,k + ), (6.28) 



where p% is the nuclear baryon density and N(k±, k + ) is the occupation number for a nucleon 
of momentum (k±, k + ). Recall that the variable k + is defined in Eq. ( 4.25| ). Since the integral 



gives the total plus-momentum carried by nucleons, the integrand (over k + ) which multiplies 
the factor k + is the probability f{k + ) that a nucleon has momentum k + distribution. Thus: 

^f = fdk + k + f(k + ), (6.29) 

where 



^ + ) = ^T^3 / d 2 k ± N(k ± , k + ). (6.30) 

A function f(y) can be obtained by replacing k + by the dimensionless variable y using 
y = ~- with M = M — 14.71 MeV. " 

The relation to experiments is obtained by recalling that the nuclear structure function 
F 2 a can be obtained from the light front distribution function f(y) (which gives the proba- 
bility for a nucleon to have a plus momentum fraction y) and the nucleon structure function 
F 2N using the relation J531: 



F 2A (x) 
A 



dyf(y)F 2N (x/y), (6.31) 



where y is A times the fraction of the nuclear plus-momentum carried by the nucleon, and 
x is the Bjorken variable computed using the nuclear mass minus the binding energy. This 
formula is the expression of the usual convolution model, with validity determined by a 
number of assumptions. Our formalism enables us to calculate the function f(y) from the 
integrand of Eq. ( |4.28| ) . 



B. Computing the total number of mesons 

We can get the total number of each kind of meson (except the u) using a sum rule. 
Consider the schematic form of the equation for the G-matrix 



G^) = V{Pr 3 ) + V (Pr)^G(Pr), (6.32) 



where is a schematic representation of the propagator of Eq. ( |4.53| ). Differentiating with 
respect to PZ yields the result: 

~dP^~ " ( AE ] ^P~ {1 + ~KE (6 - 33) 

The Moeller operators appear to the right and left of the derivative of the potential. Fur- 
thermore, in our one boson exchange approximation, the total nucleon-nucleon potential is 
the sum of the contributions due to individual bosons. Thus we may define 
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dG m (Pr 



dP~ 



1 + G 



Q , dV m {P- ) Q 
AE dPF, 1 AE 



G), 



(6.34) 



in which the label m refers to the type of meson. But The potential V m appearing in 
Eq. (|6.32| ) is simply the Fourier transform of the potential V m (q). Thus an examination of 
eqns. (6.25) and (6.24) shows that, considering the pion for example, 



AL 



d 3 qn^(q) = ^ {a,(3 

a,/3<F 



dPn 



a, (3), 



(6.35) 



Numerical evaluation of Eq. ( |6.35| ) leads to the result that N^/A = 0.05. This is smaller 
than the 18% of Friman et al |)4] because we use scalar mesons instead of intermediate A 
states to provide the bulk of the attractive force. 

The expression for the density of vector mesons involves the removal of the effects of the 
instantaneous exchange and one must use explicit light front variables. 



VII. IMPLICATIONS FOR LEPTON-NUCLEUS DEEP INELASTIC 
SCATTERING AND THE NUCLEAR DRELL-YAN PROCESS 



The values of Us, U v and N n allow us to assess the deep inelastic scattering of leptons 



from our version of the ground state of nuclear matter. Using M*(kp) = 744 MeV |y| and 
neglecting the influence of two-particle-two-hole states to approximate }'{k + ) [|],||,|56| shows 
that nucleons carry 81% (as opposed to the 65% of mean field theory [|IJ) of the nuclear plus 
momentum. This represents a vast improvement in the description of nuclear deep inelastic 
scattering. The minimum value of the ratio F2a/F 2 n-, obtained from the convolution formula 
( |6.31| ) is increased by a factor of twenty towards the data as extrapolated in Ref . |57|] . But 



this calculation provides only a lower limit of the nucleon contribution because of the neglect 



of effects of the two-particle- two hole states [58 



Turn now to the experimental information about the nuclear pionic content. The Drell- 
Yan experiment on nuclear targets |59| showed no enhancement of nuclear pions within 
an error of about 5%-10% for their heaviest target. No substantial pionic enhancement 
is found in (p,n) reactions pD| . Understanding this result is an important challenge to 
the understanding of nuclear dynamics [fPH . Here we have a good description of nuclear 
dynamics, and our 5% enhancement is consistent [32|, within errors, with the Drell-Yan 
data. 



VIII. SUMMARY AND DISCUSSION 

This paper contains a new relativistic light-front theory of nuclear matter. Light front 
quantization is used to obtain a nucleon-nucleon potential which yields phase shifts in good 
agreement with data. We use this as input in a light front many body theory. A straight- 
forward derivation leads to a theory in which the effective interaction is the light-front 
G-matrix. We obtain a good description of the binding energy, density, and incompressibil- 
ity of nuclear matter. The binding energy per nucleon is 14.7 MeV and kp = 1.37fm _1 . The 
compressibility is 180 MeV. 
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The use of a meson-nucleon Lagrangian enables us to also compute the mesonic content 
of the wave function using a consistent approach represented by Eqs. (|4.1|) , (|4.12| ) and fl6.2|) . 



The results are not in conflict with extrapolations of deep inelastic scattering and Drell-Yan 
data to nuclear matter. 

The omega and rho mesonic content are still to be evaluated. We believe that the 
possible nuclear enhancement of vector mesons is a promising avenue for future theoretical 
and experimental research. 
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APPENDIX A: NOTATION, CONVENTIONS, AND USEFUL RELATIONS 

This is patterned after the review of Harindranath J7| The light-front variables are defined 
by 

x+ = x° + x 3 , x~ = x° - x 3 , (Al) 
so the four- vector x^ is denoted 

x» = (x + ,x~,x ± ). (A2) 
With this notation the scalar product is denoted by 

x ■ y = ^x + y~ + ^x~y + - x 1 ■ y 1 . (A3) 



The metric tensor g^ v with jj, = (+, — , 1,2) is obtained from the usual one by using (AT) 



(i.e. g 0tJ- = g 0fJ- + g 3fl ). Then g H = g~ + = 2,g l i = —1, with the other elements vanishing. 
The term g^ v is obtained from the condition that g a ^g^ = S ai . Its elements are the same 
as those of g^ u except for g_ + = g + _ = 1/2. Thus 

x- = ^x + , x + = ^x~, (A4) 
and the partial derivatives are similarly given by 

d + = 29_ = 2— 8- = 28+ = 2-—. (A5) 
ox~ ox + 



The Bjorken and Drell |3g] convention for gamma matrices is used and 

7 ± =7°±7 3 . (A6) 

the relations 
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7 7 = 0, 7 7 7 = 4 7 , 7 7 7 = 4 7 (A7) 

can be used to simplify various computations. 

The Hermitian projection operators A± are given by 

A± = |7V = ^7V = ^±« 3 ), (A8) 
and obey the following relations 

(A ± ) 2 = A±, 7 X A ± , = A±t\ (A9) 

7 °A± = A t7 ° a 1 ^ = A T a x , (A10) 

7 5 A± = A ±7 5 7 T = 2A ± 7° = 7 T A T , (All) 

7'At = \i ± (A12) 



a j j i A + = -e ij 7 + 7 5 . (A13) 
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TABLE I. Potential parameters and predictions for the deuteron and low-energy np scattering. 
For the deuteron, the binding energy B^, the D-state probability Pjj, the quadrupole moment 
Qd, and the asymptotic D-state over S-state ratio D/S are given. Low-energy np scattering is 
parametrized in terms of a np and r np in 1 Sq and at and rt in 3 Si, where a denotes the scattering 
length and r the effective range. The nucleon mass is M = 938.919 MeV. 

Light-Front OBEP Thompson OBEP a Empirical 6 



Meson Parameters: 





m a (MeV) 


gl/^ [f/g] 


A Q (GeV) 


gl/^ \f/g] 


A a (GeV) 


if/g] 




138.04 


14.0 


1.2 


14.6 


1.2 


13.5 - 14.6 


V 


547.5 


3 


1.5 


5 


1.5 


< 5 


p 


769 


0.9 [6.1] 


1.85 


0.95 [6.1] 


1.3 


0.6(1) [6.6 ± 1.0] 




782 


24.5 [0.0] 


1.85 


20.0 [0.0] 


1.5 


24 ± 5 ± 7 


a 


983 


2.0723 


2.0 


3.1155 


1.5 




a 


550 


8.9602 


2.0 


8.0769 


2.0 





Deuteron 



B d (MeV) 
Pd (%) 
Qd (fm 2 ) 
D/S 



2.2245 
4.53 
0.270 c 
0.0250 



2.2247 
5.10 
0.278 c 
0.0257 



2.224575(9) 

0.2860(15) 
0.0256(4) 



Low-Energy np Scattering 



a np (fm) 
r np (fm) 
at (fm) 
n (fm) 



-23.745 
2.671 
5.494 
1.856 



-23.747 
2.664 
5.475 
1.828 



-23.748(10) 
2.75(5) 
5.424(4) 
1.759(5) 



"Potential B of Brockmann and Machleidt [|30| . 

b For more comprehensive information on the empirical data and references, see Table 4.1 



and 4.2 of Ref. 
c Meson-exchange current contributions not included. 
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FIGURES 



FIG. 1. Phase shifts 5 and mixing parameters e of neutron-proton scattering for partial 
waves with J < 2 and laboratory kinetic energies T[ a b < 300 MeV. The solid line is the 
prediction by the LF OBEP presented in Sec. Ill and the dotted line the one by Potential B 
of Brockmann and Machleidt |5(J. The open circles represent the multi-energy np analysis 
by the Nijmegen group [H3l and the solid dots are the VPI analysis SM97 §44]. 



FIG. 2. Energy per nucleon in nuclear matter, E/A (in units of MeV), as a function of 
density expressed by the Fermi momentum kp (in units of fm^ 1 ). The solid line is our pre- 
diction using light-front Brueckner theory. The dotted curve is obtained when the medium 
effect on meson propagation is omitted. The dashed line is the result from conventional 
(nonrelativistic) Brueckner theory. The box describes the area in which nuclear saturation 
is expected to occur empirically. 
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